Theory of hypothetical ferroelectric superlattices incorporating head-to-head and 

tail-to-tail 180° domain walls 
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While electrical compatibility constraints normally prevent head-to-head (HH) and tail-to-tail 
(TT) domain walls from forming in ferroelectric materials, we propose that such domain walls could 
be stabilized by intentional growth of atomic layers in which the cations are substituted from a 
neighboring column of the periodic table. In particular, we carry out predictive first-principles 
calculations of superlattices in which Sc, Nb, or other substitutional layers are inserted periodically 
into PbTi03. We confirm that this gives rise to a domain structure with the longitudinal component 
of the polarization alternating from domain to domain, and with the substitutional layers serving as 
HH and TT domain walls. We also find that a substantial transverse component of the polarization 
can also be present. 
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PACS numbers: 61.50.Ks, 77.84.Dy, 81.05.Zx 

In recent years, first-principles calculations have been 
used very successfully to explore new superlattice struc- 
tures based on ABO3 perovskites and to predict their 
physical properties Because of the vast number of 
combinations that can be assembled from the elemen- 
tal building blocks, first-principles calculations can pro- 
vide a useful way to screen for desired physical prop- 
erties and to study interesting possibilities theoreti- 
cally before they are grown experimentally. 5,6 For ex- 
ample, Sai et al. predicted that a compositional per- 
turbation that breaks inversion symmetry could in- 
crease the piezoelectric response^ and an enhancement 
of the spontaneous polarization in BaTiOs/SrTiOs su- 
perlattices was predicted by Neaton et al. under cer- 
tain conditions of epitaxial strain^ Also, the break- 
ing of tetragonal symmetry has been predicted for cer- 
tain BaTi03/SrTiC>3 superlattices^ and new phase di- 
agrams have been predicted for superlattices compose 
of Pb(Sco. 5+^^0.5-^)03 with certain layer sequences of 
composition /i£ 

In the present work, we propose a new kind of su- 
perlattice structure which might be grown artificially 
via controlled layer-by-layer epitaxy techniques such as 
molecular-beam epitaxy (MBE)1 or pulsed laser depo- 
sition (PLD)£ These contain 180° ferroelectric domain 
walls, but of an unusual kind. Normally, 180° domain 
walls in a tetragonal ferroelectric such as PbTiC>3 form 
so as to obey electrical compatibility, leading to 180° do- 
main walls running parallel to the tetragonal (polariza- 
tion) axis. This occurs because a tilting of the 180° do- 
main wall would lead to a polarization-induced bound 
charge on the domain wall, leading to a large and unfa- 
vorable electrostatic energy. Forcing the introduction of 
a "head-to-head" (HH) or a "tail-to-tail" (TT) domain 
wall perpendicular to the polarization direction would 
presumably cause it to become strongly metallic due to 
neutralization by free carriers. 

However, we propose that intentional insertion of 
chemical substitutions in certain atomic layers during 
growth ("delta doping") could provide an approximate 



cancellation of this bound charge, allowing for the fabri- 
cation of new kinds of superlattice structures containing 
180° HH and TT domain walls, as illustrated in Fig. □ 
In this PbTiC>3 example, the compensation is realized 
by substituting certain layers of Ti atoms by donor or 
acceptor atoms drawn from neighboring columns of the 
Periodic Table. If the substitutional charge density can 
be chosen to provide an approximate cancellation of the 
bound polarization charge of the HH or TT domain wall, 
then one may hope to arrive at a configuration that is 
both neutral and energetically favorable. 

In this paper, we carry out first-principles density- 
functional calculations of superlattice structures such as 
the 8-cell one shown in Fig. as a "proof of concept" 
for this idea. We find that this arrangement, with HH 
and TT domain walls at the acceptor and donor layers, 
respectively, is a robust configuration in which there is 
no gap closure (i.e., the superlattice remains insulating). 
Somewhat surprisingly, we also find that a small trans- 
verse component of the polarization develops, leading to 
an overall spontaneous polarization of the supercell in the 
[100] direction. This can be understood as arising from 
the imperfect compensation in this case. 

We first select a host material for our superlattice, 
choosing a tetragonal material and letting its c axis be 
in the superlattice growth direction. In order to attain 
perfect compensation by a donor or acceptor layer of 
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FIG. 1: Sketch of PbTiOs-based ferroelectric superlattice in 
which head-to-head and tail-to-tail 180° domain walls coin- 
cide with acceptor (Sc) and donor (Nb) substitutional layers 
respectively. 
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areal charge density ±e/a 2 , where a is the in-plane lat- 
tice constant, we would need a material whose sponta- 
neous polarization P s is ±e/2a 2 . Roughly, a ~ 4 A for 
most perovskites, in which case the matching condition is 
P s = 0.50 C/m 2 . Taking KNb0 3 , BaTi0 3 , and PbTi0 3 
in their tetragonal phases, these have P s — 0.40 C/m 2 , 
0.27 C/m 2 , and 0.75 C/m 2 , respectively. None of them is 
ideal, but we have chosen to carry out our pilot calcula- 
tions here using PbTiOa, whose polarization is somewhat 
larger than the nominal matching value. PbTiOa is a 
good test system, being a well-studied prototypical ferro- 
electric material having a simple phase behavior, remain- 
ing in the tetragonal ferroelectric phase from T = up 
to well above room temperature. Experimentally, alloy 
host materials such as Pbi-^Sr^TiOs or PbZr1_.jTi.cO3 
might be used to tune the concentration to achieve more 
ideal matching of the polarization, but we have decided 
to limit ourselves to pure compounds for this first study. 

We construct a 40-atom supercell by arranging 8 unit 
cells of bulk tetragonal PbTiC-3 in an a x a x 8c super- 
cell, with initial atomic coordinates chosen to be those 
of the relaxed bulk ferroelectric structure with polariza- 
tion pointing in the [001] direction throughout the su- 
percell. We then substitute the Ti atoms in the first 
Ti02 layer by donor atoms, and the Ti atoms in the 
fifth TiC>2 layer by acceptor atoms, as shown in Fig. ^ 
We choose Nb (5+) as the donor and Sc (3+) as the 
acceptor because these atoms are commonly found on 
the B site of other perovskite ABO3 oxides. We assume 
perfect control of the layer-by-layer composition, result- 
ing in an ideal PbNb0 3 /3xPbTi03/PbSc0 3 /3xPbTi03 
structural sequence as shown Fig. ^ 

We then relax the structure under a symmetry con- 
straint that preserves the tetragonal PAmm symmetry 
(i.e., atom displacements arc allowed only in the z direc- 
tion). The lattice constant is fixed in the in-plane direc- 
tions, corresponding to conditions of epitaxial growth on 
a bulk tetragonal PbTiOs substrate. We also fixed the 
lattice constant along the [001] direction, assuming that 
the volume relaxations resulting from the compensation 
would be small. All calculations were performed at zero 
temperature. 

The first-principles calculations were carried out us- 
ing the Vienna ab-initio simulations package (VASP)£, in 
which density-functional theory within the local-density 
approximation^ is implemented. A plane-wave basis set 
and projector augmented-wave potentials^ were used 
with a 400 eV energy cutoff for the electronic wave- 
functions. Brillouin zone integrations were performed 
on a Monkhorst-Pack mesh of dimensions 4x4x2 and 
10x10x2 for structural relaxation and density-of-states 
(DOS) calculations respectively. The ionic relaxation was 
considered to be achieved when the Hellmann-Feynman 
forces on the ions were less than 3meV/A. 

We find that the structural relaxation converges to- 
wards a structure in which the polarization reverses in 
the top half of the unit cell and the system acquires an 
M z mirror symmetry across the planes containing the 



donor and acceptor atoms, converting these layers into 
180° TT and HH domain walls as anticipated in Fig. ^ 
When we observed this apparent convergence towards the 
higher-symmetry PA/mmm structure, we carried out the 
final relaxation with this PA/mmm symmetry enforced, 
and confirmed that the resulting structure had a lower 
energy. 

In order to test whether this PA/mmm is really the 
structural ground state, we gave additional small dis- 
placements to the atoms in order to break the tetrag- 
onal symmetry, and carried out further relaxations. Sur- 
prisingly, we found that the system is unstable to the 
growth of an in-plane component of the polarization. We 
tested both initial [100] displacements and initial [110] 
displacements (preserving, in both cases, the M z symme- 
try), giving rise to Pmm2 and Amm2 structures, respec- 
tively. We find that the relaxed Pmm2 structure has an 
energy 0.065 eV lower than that of the PA/mmm struc- 
ture, while the relaxed Amm2 one is 0.057 eV lower than 
the PA/mmm. Thus, we find that the Pmm2 structure 
is our final structural ground state. 

To get a better sense of the electronic structure of the 
relaxed structures, we computed the electronic density 
of states (DOS). The results, shown in Fig. |2ta) and (b) 
for the PA/mmm and Pmm2 structures respectively, in- 
dicate that the compensation is successful and the en- 
tire superlattice is insulating. (The DOS for the Amm2 
structure looks very similar to these.) The final gap of 
1.46 eV is moderately smaller, but not drastically smaller, 
than the gap of 1.78 eV for bulk tetragonal PbTiOs. (As 
these gaps are computed using the LDA, they should be 
regarded as underestimates of the true gaps.) 

We also checked what happens if there is no compensa- 
tion between the polarization-induced bound charge and 
the donor/acceptor layer. To do this, we constructed a 
40-atom PA/mmm supercell similar to that of Fig.[5Ja), 
but without introducing the impurity layers. (Actually, 
this was done by imposing the displacement pattern of 
the ferroelectric tetragonal ground state of PbTiOs, but 
with opposite signs in the bottom and top halves of the 
supercell.) The result is shown if Fig. Efc), where it is 
clear that the structure is strongly metallic. Thus, the 
compensation by the donor and acceptor layers is essen- 
tial in order to arrive at an overall insulating structure. 

We return now to the pattern of ground-state atomic 
displacements that we calculated for the PA/mmm and 
P2mm structures of Fig. |2Ia) and (b), and interpret these 
in terms of a local polarization in each unit cell. To do 
this, we model the local polarization asi^ 

c a 

where i is the Pb-centered unit cell index, a runs over all 
atoms, and w a is a weight factor that is unity inside the 
cell, zero outside, and takes appropriate fractional values 
for atoms shared with a neighboring cellii^ For the dy- 
namical effective charge tensors Z* , we use bulk PbTiOs 
values of Z*(Pb)=3.90, Z*(Ti)=7.06, Z,T(0)=-5.83, and 
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FIG. 2: Density of states (DOS) of 40-atom PbTi0 3 su- 
percells considered in the text, (a) PA/mmm supercell of 
Fig. Q containing HH and TT domain walls compensated by 
impurity layers, (b) Same, but with in-plane component of 
polarization (Pmm2 symmetry), (c) Reference structure in 
which polarization reverses as in Fig. but no compensating 
impurity layers are present. 



Zj_(0)=— 2.56 following Ghosez et al>i£ We also arbitrar- 
ily adopt Z*(Nb)=8.06 and Z*(Sc)=6.06 based on the 
rough expectation that these atoms should resemble Ti 
except a shift by ±e of the core charge. 

The results of this analysis are shown in Fig. [3J which 
shows the local polarizations in the ground-state Pmm2 
structure of the superlattice. The local polarization in 
the [001] direction changes sign across the dopant layers, 
confirming that HH and TT domain walls are formed. 
In addition, there is a roughly uniform polarization in 
the [100] direction, although it is slightly enhanced in 
the vicinity of the HH and TT domain walls where the 
[001] polarization is crossing through zero. Thus, the en- 
tire superlattice structure is antiferroelectric as concerns 
the [001] polarizations, but is ferroelectric in the [100] 
direction. (The P z profile of the P4/mmm-constrained 
structure is not shown, but it is almost identical to that 
of the Pmm2 structure.) 

In order to understand our results, recall that PbTiC>3 
has a spontaneous polarization of 0.734e/a 2 , somewhat 
larger than the ideal value of 0.5e/a 2 that would give 
perfect compensation of the donor and acceptor layers. 
If each domain were to develop its full bulk sponta- 
neous polarization, there would be substantial charges 
(±0.47e/a 2 ) at the donor layers due to overcompensa- 
tion by the domain walls. The resulting depolarizing field 
would reduce P z to a value much closer to, but still larger 
than, the reference value of 0.5e/et 2 . This is precisely 
what we see in Fig. [3J in the middle of each domain, P z 
attains a value just slightly larger than 0.5e/a 2 . 

The appearance of the polarization in the in-plane di- 
rection can be understood based on the concept of polar- 




-0.8 - 

Pb centered cell index of superlattice 



FIG. 3: Local polarization profile of the relaxed Pmm2 struc- 
ture. The [001] component changes sign at the dopant layers, 
while the [100] component is roughly uniform throughout the 
supercell, giving rise to an overall polarization in the [100] 
direction. Dotted lines indicate the "ideal" value 0.5e/a 2 ; 
dashed lines represent the spontaneous polarization of bulk 
tetragonal PbTi0 3 . 

ization rotation. Many previous worksi^ have provided 
support for a picture in which the crystal is regarded as 
having a preferred magnitude of the polarization, with a 
rather weak crystalline anisotropy determining its direc- 
tion. If the polarization is prevented from attaining its 
preferred bulk value in the [001] direction, it can still re- 
main close to this preferred magnitude by tilting off this 
axis. Indeed, the resulting magnitude of the polarization 
in the middle layers of the slab, 0.57e/a 2 , is only slightly 
lower than the preferred bulk value. Moreover, this pic- 
ture naturally helps explain why the in-plane component 
of the polarization grows near the HH and TT domain 
walls, where P z is in the process of crossing through zero. 

Fig. 21 gives further evidence of the physical picture 
outlined above. Here, we have plotted the site-projected 
DOS for Ti atoms sitting on three different TiC>2 lay- 
ers of the Pmm2 supercell. Specifically, Ti(Nb) is a Ti 
atom sitting next to the donor-compensated (NB) plane; 
Ti(Sc) sits next to the acceptor-compensated (Sc) plane; 
and Ti(Mid) is a Ti atom midway between the two im- 
purity layers. The zero of energy is at the valence-band 
maximum. The vertical dashed line provides a reference 
that helps clarify that the spectral weight of the conduc- 
tion band, and in particular its low-energy edge, shifts to 
higher energies as one progresses from the Ti(Sc) to the 
Ti(Nb) layer. The smallness of these shifts confirms that 
the compensation is quite good, and the small shift that 
does occur can be attributed to the leftover depolariza- 
tion field in the [001] direction. Furthermore, the sign of 
this effect is consistent with our previous conjecture that 
the domain walls are slightly under-compensated. 

It is interesting to ask what will happen for super- 
lattices of longer and longer period in the c direction, 
for a given material system. Unless the compensation is 
perfect, the depolarizing fields resulting from the imper- 
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FIG. 4: DOS projected onto three different Ti atoms in the 
Pmm2 supercell. Ti(Sc) denotes the Ti atom in the Ti02 
iayer ciosest to the Sc02 iayer; Ti(Nb) denotes the Ti atom 
in the Ti02 iayer ciosest to the Nb02 layer; and Ti(Mid) 
denotes a Ti atom in the middle of the supercell, between the 
two dopant layers. 

feet compensation should result in an electrostatic po- 
tential having a sawtooth dependence on z. This should 
be sufficient to close the gap, and drive the superlatticc 
structure metallic, when the amplitude of this sawtooth 
potential exceeds the band gap. We can estimate the 
critical superlattice period for the PbTiOs/Sc/Nb sys- 
tem by recalling that we estimated a net bare charge of 
2P S - e/a 2 ~ 0.47e/a 2 or about 0.51 C/m 2 . Taking the 
theoretical value of the bulk band gap of E g = 1.8 eV and 
a bulk dielectric constant of about 110, we estimate that 
the depolarization field would close the band gap when 
the distance L between impurity layers exceeds a critical 
value of about 75 A or about 19 PbTiC>3 layers. 

If it is desired to increase the superlattice period and 
still keep an insulating structure, then the precision of 
the compensation should be improved. Two approaches 
suggest themselves. First, one could make use of mixed 
dopant layers of composition Tii_ x Sc r or Tii-^Nb^, 
with x tuned to obtain the correct degree of compensa- 
tion for a given bulk material. Second, the bulk material 
itself may be varied in order to adjust its spontaneous 
polarization closer to the ideal value of 0.5 e/a 2 . For ex- 
ample, replacing PbTi0 3 by PbZr^Tii^Oa (PZT) with 



x ~ 0.2 may accomplish this^ Using such strategies, it 
may be possible to construct long-period insulating su- 
perlattices, perhaps even with periods large enough to be 
of use for optical grating applications. 

In summary, we have calculated the properties of a 
particular realization of a hypothetical class of ferroelec- 
tric superlattice structures containing head-to-head and 
tail-to-tail 180° domain walls compensated and stabilized 
by the intentional insertion of atomic charged-impurity 
layers. In the system we chose for study, in which Sc 
and Nb substitutional layers provide the compensation 
in PbTiC>3, the spontaneous polarization of the bulk ma- 
terial is larger than the ideal value; consequently, the su- 
perlattice develops an in-plane ferroelectric polarization, 
while also developing the expected antiferroelectric do- 
main structure in the growth direction. This appears to 
be a novel way of forcing polarization rotation into the in- 
plane orientation in PbTiC>3, and it would be interesting 
to explore how the superlattice properties (spontaneous 
polarization, Curie temperature, etc.) could be tuned by 
the growth conditions (e.g., superlattice period). More 
generally, it would be interesting to investigate other host 
perovskite ferroelectric materials, other kinds of compen- 
sating layers (including, perhaps, dopants from column II 
or VI of the periodic table), and alloy compositions of the 
bulk and/or dopant layers. With a more suitable match- 
ing of the compensation, the driving force for polarization 
rotation could be reduced, and a true PA/mmm super- 
lattice structure could presumably be stabilized. In fact, 
it also seems possible that a normally rhombohedral bulk 
ferroelectric like BaTiOa whose spontaneous polarization 
is smaller than the ideal value of 0.5 e/a 2 could have its 
polarization rotated in the opposite way, into alignment 
with the z-axis, since in this case the compensating layers 
would enhance the polarization in that direction. Finally, 
much longer-period superlattices might be obtained by 
very accurate matching of the compensation. There are 
thus many open avenues for future theoretical and exper- 
imental studies. 
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